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ABSTRACT 


Early-type galaxies (ETGs) are known to be enhanced in a elements, in accordance with their old ages and short 
formation timescales. In this contribution we aim to resolve the enrichment histories of ETGs. This means we study the 
abundance of Fe ([Fe/H]) and the a-element groups ([a/Fe]) separately for stars older than 9.5 Gyr ([Fe/H]o, [a/Fe]o) 
and for stars between 1.5 and 9.5 Gyr ([Fe/H]i, [a/Fe]i). Through extensive simulation we show that we can indeed 
recover the enrichment history per galaxy. We then analyze a spectroscopic sample of 2286 early-type galaxies from the 
SDSS selected to be ETGs. We separate out those galaxies for which the abundance of iron in stars grows throughout 
the lifetime of the galaxy, i.e. in which [Fe/H]o < [Fe/Hji. We call those consistent with self-enrichment, while the 
others must have experienced some mergers or signiheant gas accretion. We conhrm earlier work where the [Fe/H] and 
[a/Fe] parameters are correlated with the mass and velocity dispersion of ETGs. We emphasize that the strongest 
relation is between [a/Fe] and age. This relation falls Into two regimes, one with a steep slope for old galaxies and 
one with a shallow slope for younger ETGs. The vast majority of ETGs in our sample do not show the ‘knee’ in the 
plot of [Fe/H] vs. [a/Fe] commonly observed in local group galaxies. This implies that for the vast majority of ETGs, 
the stars younger than 9.5 Gyrs are likely to have been accreted or formed from accreted gas. The properties of the 
intermediate-age stars in accretion-dominated ETGs indicate that mass growth through late (minor) mergers in ETGs 
is dominated by galaxies with low [Fe/H] and low [a/Fe[. The method of reconstructing the stellar enrichment histories 
of ETGs introduced in this paper promises to constrain the star formation and mass assembly histories of large samples 
of galaxies in a unique way. 


1. Introduction 

The stars in massive early-type galaxies (ETGs) have had 
to form early in the history of the universe for two reasons: 
1) Very few of the se galaxies have (significant) young stellar 
populations (e-g-, Vi et al.|2005 Trager fc Somerville|2009 


Thomas et al. 2010p and 2) the ratio of a elements (O, 

Mg, etc.) to F'e in the atmospheres of these stars is higher 
than solar (e.g., Peterson 1976 Worthey et al.|19'^ Milone 
et al.|2000|), indicating a short timescale for the enrichment 
of the ISM out of which they formed. The ratio [g/Fe] is a 


powerful estimator of the duration of star formation (jTins- 


1979 Matteucci & Greggio 


19861. Study of the mass 


dependence of the ages and enhancement ratios of ETGs 
has led to the classical picture in which the star f ormation 
timescales are shorter with increasing galaxy mass (Thomas 

eral]|200^ . 


Direct studies of massive galaxies at high redshift sup¬ 
port this picture in which at least some ETGs must have 
"formed the bulk of their mass only a few million years af¬ 
ter the big bang" (Nayyeri et al. 20141. On the other hand 
there is evidence that massive ET Gs cannot be described by 
a stellar populatio n with a single age and abu ndance pat- 
tern (recently, e.g., Lonoce et al.|2014|. Indeed, Arimoto & 


Yoshii (19861 emphasize that this is true in terms not only 
of an age spread, but also of an abundance spread: "The 
present analysis shows that such metal-poor stars in the 
RGB evolution have a strong influence on the integrated 
colors of the galaxy." 

Recent observational and theoretical work has focused 
on the question of when and how the late mass growth of 
early-type galaxies takes place. Two opposite possibilities 
for such mass growth exist: 1) in situ star formation from 
gas that was accreted in the initial phase of the galaxies for¬ 
mation and 2) mass growth from later major or minor dry 
merging. A third possibility with less clear observational 
signatures is 3) stellar mass growth through star formation 
from gas that was either accreted during a merger (wet 
mergers) or in a cold mode from the cosmic web. 

Different observational approaches to these questions 
exist. An important piece of the puzzle is the observation 
that the progenitors of ETGs at redshifts above z~l are 


compact, dense objects dTrujillo et al. 

2006 van Dokkum 

et al. 2008, van der Wei et al.||2008||. 

JNaab et al. ( 

2009), 

Hopkins et al. 

(2009p, and van Dokkum et al. (|2010a 

1 show 


ing the dense galaxies observed at high redshift, although 


van der Wei et al. (20111 make the point that the ma- 
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jority of these distant, dense galaxies should be thought 
of as disk-like, introducing some tension with the mostly 
dispersion-dominated cores of today’s early-type galaxies. 
Overall, the data have led to the notion that (minor) merg¬ 
ing over a Hubble time will simultaneously contribute to 


dGlemens et al.||2009l Naab et al.||2009 van Dokkum et al. 

2010a| User et al.||2012 Greene et al.||2012| 2013||. 

Simulati 
tant at larg 

Dns predict that stellar accretion is most impor- 

3 radii (Lackner et al. 2012 Navarro-Gonzalez 

et al. 

2013 

Hirschmann et al.||2014||. In contrast to earlier 

notions, sim 
is needed to 
On the othe 
is available 

ulations even show that no major merger at all 

form massive ellipticals (Bournaud et al.|2007l. 

r handj direct measurement of major merging 

n the literature (e.g.. Bell et al.||2006 de Ravel 

et al. 

2009 

Robaina et al. 2010 |Robotham et al. 2014p. 


dergone one major merger after redshift one (Keenan et al. 


20141 and that the higher the mass of a galaxy, the more 


it grows through major merging. The evolution in the mass 
function can be combined with the expected growth due to 
star formation to yield a d irect measure of the mass growth 
due to mergers in ETGs (Walcher et al. 20081 Ownsworth 
et al. 20141. This leads to estimates of roughly half of the 


present-day total stellar mass of ETGs being formed in 
other galaxies prior to their incorporation into the ETG 
over cosmic time. However, merge rs also contribute gas. 
For example Kormendy et al. (20091 discuss how structural 
properties can be used to distinguish between dry merg¬ 
ers that would contribute only stars and wet mergers that 
would in turn contribute gas, thus leading to late star for¬ 
mation and rebuilding of cusps. 


The present paper is dedicated to measuring the rel¬ 
ative importance of the mass-growth processes from their 
chemical signature. In particular, other observational ap¬ 
proaches cannot hope to distinguish between in situ forma¬ 
tion of stars from gas left over from any putative "initial 
accretion", on the one hand, and gas subsequently accreted 
either in a cold mode or through minor mergers, on the 
other. This distinction is much easier to study from chem¬ 
ical signatures as shown below. 

In this contribution we define the "stellar enrichment 
histories" (SEH) as the four-dimensional space occupied by 
the stellar populations in each galaxy and spanned by the 
parameters age (i.e., time elapsed since formation of those 
populations), [Fe/H], [a/Fe] and the mass (or luminosity) 
contribution of these stars. It needs to be kept in mind 
that we can only address the stellar population within one 
spatial resolution element (one fiber in our case, the case 
of SDSS). Because of measurement uncertainties we very 
strongly bin this space into just two age bins: the old age 
bin contains all stars older than 9.5Gyr, the intermediate 
age bin contains all stars between 1.5 and 9.5 Gyr of age. 
Within each age bin and for each galaxy we can compute 
the luminosity-averaged [Fe/H] and [a/Fe] as well as the 
fractional contribution of the old and intermediate popula¬ 
tions to the total luminosity of each galaxy. This principle is 
illustrated in the upper panel of Figured For the remainder 
of this paper we denote the mean [a/Fte] and [Fe/H] of the 
old stars with the subscript ’o’ (i.e. we write [Fe/H[o and 
[Q;/Fe[o), whereas we denote the mean [a/Fe] and [Fe/H] of 
the intermediate age stars with the subscript ’i’ (i.e. we 
write [Fe/H[i and [a/Fe[i). The "star formation history" 
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(SFH) would simply be the marginalization of the SEH onto 
the 2-dimensional space of age and luminosity contribution. 


The choice of 9.5 Gyr for separating old and intermedi¬ 
ate stellar populations in ETGs is of course to some level 
arbitrary. It can, however, be justified on several grounds 
as providing the widest age bin for old stars that still makes 
se nse physical l y. On the one hand simulations such as those 
of Naab et al. (20091 show the typical mass assembly histo- 
ries of ETGs in a hierarchical universe. For many of these 
galaxies there is a change of regime at around 10 Gyr ago, 
where the mass build-up changes from being dominated by 
in-situ star formation to accretion of other, smaller galaxies 
and clumps. Thus this is an age where we expect to maxi¬ 
mize the signal from any difference between old and inter¬ 
mediate age population. On the other hand, those alpha- 
enhanced populations that we can study in detail in the 
Milky Way are all older than II Gyr. This is true for the 


thick disk (Fuhrmann 20111 and for the bulge (Bensby et al. 


20131. As shown later in the present paper (Figure 12, lower 


right panel) and in a paper in preparation (Walcher et al., 
in prep.), this effect is seen in the luminosity-weighted av¬ 
erage properties at a lookback time of around 9.5 Gyr as 
well in the sense that the rate of change of [a/Fe] with 
time shows a significant change of slope. These 2 or 3 Gyr 
after the onset of star formation are thus different in the 
chemical signatures, which helps for the technical as well as 
scientific arguments of the present paper. Finally, the cos¬ 
mic star formation density peaks at a lookback time of 10 
Gyr, a s shown in the recent review by [Madau fc: Dickinson] 
(2014| their Figure 9). Specifically for ETGs, this is also 
the time when qu iescent galaxies start ed appearing in large 
number (see e.g., Arnouts et al. ][2007l , their Figure 13). 

The lower panel of Figure illustrates basic ex¬ 
pectations for the locus the old and intermediate stel¬ 
lar population would occupy in the projection onto the 
[a/Fe] vs. [Fe/H] plane. As is well known the ’knee’ in 
this diagram carries information about the star formation 
efficiency in the galaxy, where higher star formation effi¬ 
ciency leads to higher overall metallicity before the onset 
of SNIa enrichment, and thus to a knee that is located at 
higher [Fe/H] values. The low-[Fe/H] plateau in [a/Fe] can 
be modulated by the slope of the upper end of the stellar 
initial mass function, but the data we discuss here do not 
allow any inference concerning this, thus we neglect this 
question for the remainder of the paper. 


The existence of the knee in the lower panel of Figure [^ 
for the Milky Way stars is tied to some causal connection 
at each time step during the star formation history between 
the existing stars and the next generation of stars. Although 
some accretion is needed, if only to solve the g-dwarf ques¬ 
tion, we generally call this a scenario that is consistent with 
self-enrichment (SLF). The SLF scenario would predict in 
our specific case of coarsely resolved enrichment histories for 
ETGs that the [Fe/H] of the intermediate-age population 
would always be higher than for the old population and the 
[a/Fe] of the intermediate-age population would always be 
lower than that of the old population. On the other hand, as 
mentioned above we know that mergers contribute a signifi¬ 
cant percentage of the stars in massive ETGs. The merging 
stars bear no connection to the pre-existing stars in the 
ETG and thus could have any [Fe/H] and [a/Fe] proper¬ 
ties. The shorthand for this scenario throughout this paper 
is AGG for accretion. 
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Lookback time [Gyr] 



[Fe/H] 


Fig. 1. Schematic representation of the new parameters intro¬ 
duced in this paper, namely [a/Fe]o and [Fe/H]o. These repre¬ 
sent the average [a/Fe] and [Fe/H] of stars older than 9.5 Gyr 
on a per galaxy basis. The upper panel shows a schematic star 
formation history, which separates into these two stellar pop¬ 
ulations. The lower panels depicts on the standard [a/Fe] vs. 
[Fe/H] diagram familiar for the Milky Way what is expected for 
the old and intermediate age stars respectively. 


For example, Pipino et al.|(2009|) show that semi-analytic 
models cannot reproduce the [a/Fe[-mass relations yet, ow¬ 
ing to inadequate treatment of late mergers, thus clearly 
implying that we do not have sufficient theoretical under¬ 
standing of low-mass galaxy evoluti on in simulations. More 
generally, Hirschmann et al. (20131 conclud that, while it 
is clear that accretion, minor merging, major merging, star 
formation efficiency, and winds must all work together to 
form those galaxies we see today, their exact interplay is 
still not sufficiently understood. Simulations thus do not 
produce a satisfying population of galaxies overall. Obser- 
vationally, and in the direct context of this contribution, 
this is illustrated by the recent work of |Renzini fc Andreon] 
(20141, who concluded from data on the enrichment of the 
intracluster medium "that even the most massive galaxies 
must have lost a major fraction of the metals they have 
produced". 

An uncertainty that much more directly affects our 
distinction in SLF and ACC ETGs has been spelled out 


in the works of Arimoto & Yoshii (19861 and Vazdekis 


et al. (19961. Indeed, it turns out that metallicity may 


not monotonously increase with age for the most massive 
ETGs even in absence of accretion events. After a peak in 
metallicity at a certain age, this metallicity starts declin¬ 
ing again because of the late mass loss by the abundant 
low-metallicity stars st ill present in the gala xy. Compar¬ 
ison with Figure 11 of Vazdekis et al. (19961 reveals that 
the peak in metallicity may be reached very rapidly (within 
2 or 3 Gyrs after formation of the galaxy). However, the de¬ 
cline in metallicity remains subtle and is of order 0.2 dex 
in [Z] until redshift zero. We come back to this value in 
Section 16.31 

The present contribution introduces extensions to the 
method of spectral fitting. We therefore present our al¬ 
gorithm in detail in Section and then spend significant 
space on understanding any possible method-intrinsic ef¬ 
fects, such as degeneracies, in Section Finally we also 
verify in detail that the method is applicable to the data 
at hand in Section Readers wishing to read about the 
scientific results first can jump directly to Section 


The mass-metallicity relation (Tremonti et al. 2004 


Gallazzi et al. 20051 leads to the expectation that mi¬ 


nor, i.e. low mass mergers would contribute stars with 
sub-solar [Fe/H[. The [a/Fe] of these stars is less certain. 
However, Sansom & Northeast (20081 show that the high 
[of/Fe] abundances seen in massive quiescent galaxies do 
not continue to the lower luminosity quiescent galaxies, as 
those have generally lower [a/Fe]. Also, lower mass galax¬ 
ies tend to have prolonged star formation histories, again 
leading to an expectation of lower [a/Fe] in general. We 
caution that we currently do not know what the [a/Fe] of 
intermediate-redshift minor mergers could be. Neverthe¬ 
less, we are led to an expectation that the intermediate- 
age stars in ETGs would on average have solar [a/Fe] and 
[Fe/H] either slightly above or even significantly be low the 
solar value (the blue ellipse in Figure [B. Indeed, Goelho 


et al. (20091 found that for the galaxy M32 the metal-poor 


component of their stellar population fit corresponds to a 
younger population than the metal-rich one. 


Although simple, this picture can be used to provide a 
framework for interpreting the results in this paper. Gom- 
plications arise, however, that need to be acknowledged. 


2. Sample 

For this paper we are interested in a significant sample 
of early-type galaxies (ETGs) with high S/N spectra. 
With ETG we mean a "spectral ETG", i.e. essentially a 
galaxy with a spectrum dominated by old stars. Such a 
sample is easil y selected from the M PA/JHU catalog^ 
for SDSS DR7 (Abazajian et al. 20091, containing generm 
information about the galaxy’s spectrum (S/N, stellar 
velocity dispersion), emission lines and absorption indices, 
as well as photometric information. ETGs in the SDSS 
have been studied in many contributions to the literature 

(e-g-, ^ " 


Bernardi et al. 2003 
these 


Zhu et al. 

ttat 


20 lO]). The specific 
hey provide 


reason to revisit these spectra is that tney provide a 
large sample of readily available high quality spectra with 
which to test the power of the new models and the new 
methodology we apply here. In particular we adopt the 
following criteria to select our sample: 


— Signal to noise larger than 40 
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30 40 50 60 70 80 90 100 

Signal-to-noise 

Fig. 2. S/N distribution of our sample galaxies. 


— r-band concentration index larger or equal to 2.8 

— "unclassifiable", i.e. with S/N <3 in all of the BPT di¬ 
agram emission lines, thus excluding AGN 

— Velocity dispersion between 40 and 375 km/s 

— Eliminate duplicate observations, keeping the one with 
higher S/N 


We further estimate stellar masses from absorption in¬ 
dices following the exact same method as [Gallazzi et al.| 
(|2005 1. The sample is thus further restricted to those galax- 
ies for which a spectroscopic stellar mass could be obtained 
(i.e. for which the set of five absorption features, r — i fiber 
color and z-band model magnitude are available). We use 
the 16, 50, 84 percentiles of the probability distribution 
function for the lower Icr errorbar, the actual value and the 
upper 1(7 error bar, respectively. These masses are normal¬ 
ized to the z-band modelmag and a Ghabrier IMF. 

These criteria give us a sample of 2286 pure early-type 
galaxies. This sample has been S/N selected and it would 
thus be a daunting task to follow through all selection cri¬ 
teria to infer how representative the sample is for all early- 
type galaxies. We thus refrain from such an exercise. The 
S/N distribution of our sample is shown in Figure]^ 


3. Stellar population models and fitting procedure 

Stellar population models exploring the effects of varying 
element abundance ratios on specific line indices have been 
drivi ng our understan ding of ETG formation and evolution 
(e.g., Lee et al. 20091. In part icular, studying the depen- 
danci^of s ingle elements (e.g.. Smith et al.|2009 Worthey 
et al.|[2M4| is interesting in its own right, but is currently 
not possible for the goals of this paper for two reasons: 1) 
We want to study SEHs in a resolved way and it would 
be difficult to control physical degeneracies. 2) We need 
spectral models with a restricted set of free parameters to 
keep our algorithms mathemati cally under con t rol. T aking 
e.g. the recent contribution by Worthey et al. (20141 as a 
comparison and referring back to Figure we can high¬ 
light one important point of our paper. These authors con¬ 
sider abundance distribution functions (ADFs) within each 
ETG to address the similar questions to the present con¬ 
tribution, namely how abundance ratios are spread within 
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ETGs. However , due t o the nature of the fitting algorithm 
Worthey et al. (2014) have to parameterize these ADFs, 
thus leading to a prior assumption on the SEHs possible 
in the galaxies they analyze. Because knowledge of the real 
SEHs of ETGs is so scarce we do not believe that putting 
such a prior is warranted yet. In support of our approach we 
also highlight that the use of [a/Fe] as a single parameter 
involving a set of elements is also justified by the findings 


of Johansson et al. (20121, who find that the correlations of 


[0/Fe|, [Mg/Fej, [G/Fe] with velocity dispersion are identi¬ 
cal for ETGs. In general, models capable of predicting the 
effects of Qf-elements on the full optical spectrum of stel- 


et al. (20071), 

K^oleva et al. (2008c|), IWalcher et al. (|2009D, 

Percival et al. 

(2009D, Vazdexis et al. 1(20111), Gonroy &: van 

Dokkum (2011: 

), Prugniel & Koleva (2012p, Vazdekis et al. 


We use an updated versi on of the differential st ellar pop¬ 
ulation models presented in Walcher et al.| (2009 hereafter 
W09) and we refer to that paper for a full description. In 
brief, to compute the models we combine the simple stellar 
population (SSP) spectra from f ully th eoretical stellar pop¬ 
ulation models by |Goelho et al. ( 2007| (based on isochrones 
from V^iissetal. 1 12007| ) with semi-empirical i nodels by 
Vazdekis et al.|p010p and Falcon-Barroso et al. (12011 1. Of 
the Vazdekis et al.| ( |2010 ) models we only use those with 
solar metallicity and abundance ratios. The model spectra 
for non-solar [Fe/H] and [a/Fe] values are then obtained 
using the G07 models of the same age as a prediction for 
the effects of varying [Fe/H] and [a/Fe[. We show in WOO 
that these models yield reliable values for the ages, the 
[Fe/H] and the [a/Fe] values of Milky Way globular clus¬ 
ters within the range covered by our models. 

During the course of the project we realized that the 


parameter coverage of the Goelho et al. (2007) models 


(age: 3-12 Gyr, [Fe/H[: -0.5, -0.25, 0.0, 0.2 and [a/Fe[: 
0.0, 0.4) lead to "edge effects", in the sense that when fit¬ 
ting a stellar population older than 12 Gyr and with solar 
abundance ratios, age-metallicity degeneracy will drive the 
best-fitting model to have super-solar [Fe/H] and [a/Fe[. 
We have therefore decided to extrapolate the grid of the¬ 
oretical models to an age range of 2-13 Gyr, [Fe/H] val¬ 
ues of —0.5, —0.25, 0.0, 0.2 and [a/Fe] values of —0.2, 0.0, 
0.2, 0.4. The theoretical spectra are interpolated linearly in 
a multi-dimensional space with the axes: pixel, log(Flux), 
age, [Fe/H] and [a/Fe] before being applied to the semi- 
empirical models. 

Metallicity is defined as a combination of [Fe/H] and 
[a/Fe], ex plicitly for ou r models [Z] = [Fe/H] + 0.75 * 
[a/Fe] (see Goelho et al.||2007 ). 

Although we select galaxies for our analysis that have 
all signs of lacking young stellar populations, we cannot 
rule out any contribution. In the literature, rejuvenation of 
ETGs is still a controversial subject with findings related 
to sample definition, analysis methods and uncerta inties re- 
lated to late, hot phases of stellar evolution ([Fe/H] [Yi et al. 


20051 


et al. 


temp 


Sanchez-Blazquez et al.|[2009| |Ocvirk||2010| [Thomas 


20101. W e therefore extend our template set by five 

([2010 ) with younger ages 

hT 


ates from Vazdekis et al. 
(1, 0.63, 0.32, 0.1, 0.06 Gyr). ' 
[a/Fe] 


lese do not have varying 
ratios. We include these templates into the fitting 
in all cases and use them to verify that indeed contributions 
by younger stellar populations are negligible in the data we 
analyze. 
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Full spect rum fitting has survived a number of stringent 
tests ([Fe/H] Cid Fernandes et al. 2005, Koleva et al. 2008a I 
and is becoming increasingly popular as a method to derive 
stellar population properties of gala xies ([Fe/H] Heavens 
appellari fc Emsellem| [2004 |Cid h'ernandes 


et al. 


et al.| 2013 to cite but a few). For analysis of the spectra 
using full spectrum fitting we use the software paradise. 
The underlying algorithm is fully described in Appendix 
[A| and the same software w as used in our previous papers 
( Walcher et al.|[^006 2009). In one sentence, paradise fits 
model template spectra to observed spectra taking into ac¬ 
count the uncertainties on each spectral pixel, solving for 
the kinematics of the stellar population and allowing very 
flexible masking of wavelength regions due to quality con¬ 


cerns. 


4. Analysis of mock galaxy spectra 

The goal of this paper is to shed new light on the SEHs of 
early-type galaxies to understand their formation history. 
Before we can do so, we must verify that our method is in 
principle capable of providing us with SEHs. In the present 
section we describe the creation of mock galaxy spectra 
and their analysis. As we use the same models for the mock 
spectra as we then use to analyze these, we can only test the 
accuracy and intrinsic degeneracies of our method under the 
assumption that the model is a perfect representation of the 
data. We performed an independent test of the accuracy 
of our models as com pared to globular cluster spectra in 
Walcher et al. (2009). In the remainder of the paper we 
report what we find when fitting the wavelength range 4828 
to 5364A. This wavelength range has been shown in W09 
to most accurately recover age, [Fe/H] and [alpha/Fe] for 
bulge globular clusters. 

The main conclusion of this section is that we can re¬ 
cover the luminosity-weighted average properties of ETGs 
and the element abundance patterns of their old stars with 
good accuracy and precision. Recovering the properties of 
intermediate-age stars is a more challenging task. 


4.1. Mock galaxy spectra 

As motivated in the introduction, we analyze the SEHs of 
our galaxies in three bins of age r, namely r < 1.5 Gyr 
(called young age bin from now on), 1.5 Gyr < r < 9.5 
Gyr (called intermediate-age bin, subscript i), r > 9.5 Gyr 
(called old age bin, subscript o). We do not attempt a better 
time resolution for the reasons given below and because ex¬ 
perience with full spectrum fitting shows that three bins in 
age are a reasonable assumption for the best age resolution 
one can hope for in the star formation history (compare e.g. 
Gid Fernandes et al.|2005 Tojeiro et al. 2007| ). We therefore 
construct our mock galaxies by separating our SSPs in the 
three age bins above and assuming a constant SFR within 
each of these bins. We then go on to neglect the youngest 
age bin, as any young stars (or hot stars, to be more pre¬ 
cise), if present, represent a very small fraction of the light 
for the larg e majority of early-type galaxies (90% Thomas 

eTZpoM. 

We then compute a grid of mock galaxies over the fol¬ 
lowing parameter ranges: [Fe/H]i and [Fe/H]o vary between 
-0.5 and 0.2, [a/Fe], and [a/Fejo very between -0.2 and 0.4 
and the mass contribution of the intermediate stars to the 
total stellar mass may vary between 0 and 100 %. We note 








Fig. 3. Simulations of parameter recovery for a randomized 
coverage of parameter space and a total of 3000 mock galaxies. 


that we thus do not assume an exponentially falling star 
formation history at any point. On the other hand, for this 
exercise we divide the SFH of our mock galaxies in two age 
bins with constant star formation rate. Our grid of mod¬ 
els is probably more varied than real early-type galaxies, 
in particular concerning the contribution of intermediate- 
age stars, and thus any conclusions we derive from here 
is applicable to a more restricted set of SFHs as well. All 
model spectra are convolved with a velocity dispersion of 
150 km/s. We verified that our results do not depend on 
velocity dispersion, as long as this parameter is correctly 
determined. One of the advantages of full spectrum fitting 
is that the velocity dispersion is determined simultaneously 
with the stellar population analysis, thus avoiding some 
of the potential pitfalls with the more classical Lick index 
method. Guided by the S/N distribution of our sample (see 
Figure]^, we perturb the spectra assuming Gaussian un¬ 
certainties to the lowest S/N ratio of the sample, i.e. 40. 


4.2. Recovery of light-weighted mean quantities 

We now use pcuradise to recover the following quantities: 

— V-band light-weighted mean parameters age, [Fe/H], 
[of/Fe], Z 

— V-band light-weighted mean parameters for old compo¬ 
nent age^, [Fe/H[o, [a/Fe[o, Zo 

— V-band light-weighted mean parameters for intermedi¬ 
ate component age,, [Fe/H[i, [a/Fe[i, Zi 

The mathematical equation defining these values as im¬ 
plemented by paradise is 


(A) = 


AT 

l-ii=0 ' 


rLYX,, 




( 1 ) 


where are the weights for N templates with properties 
Xi as well as V-band luminosities . This equation can 
also be used to define the same quantities over a restricted 
range of SSPs that are part of the SEH (e.g. those repre¬ 
senting only old or only intermediate-age stars) by suitably 
restricting the indices over which the sum is computed. 
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Fig. 5. Verification of the uncertainties derived by the fitting 
routine through bootstrap. Each A is defined as the recovered 
minus the true value, but this time divided by the size of the 
error bar (denoted ct) as output by the fitting routine. For per¬ 
fectly Gaussian errors, the distribution should coincide with a 
Gaussian of width one and normalized to the same total area 
(dotted line). 
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Figure 3), Gid Fernandes et 

al. ( 

2005 , their Figure 6), 



and 

Sanchez-Blazquez et al. 

(|201 

their Figure 7) shows 


Fig. 4. Degeneracies in the overall mean light-weighted values 
as difference plots. Each A is defined as the recovered minus the 
true value. The dominating degeneracy occurs between age and 
[Fe/H], while age and [a/Fe] show only weak degeneracy. 


Our results for the overall light-weighted quantities, 
i.e. when the index i in Eq. [^runs over all templates, are 
summarized in F igure [3| 

Kolevaet al. (2008b) and|Sanchez-Blazquez et al. (2011) 
study the inhuence of a metallicity - velocity dispersion de¬ 


generacy, implying that when both parameters are fitted 
simultaneously, the velocity dispersion may be underesti¬ 
mated simultaneously with the metallicity or [Fe/H]. We 
find the opposite effect in our simulations, i.e. for lower 
values of [Fe/H] velocity dispersion will be over-estimated 
slightly (for our fiducial value of 150 km/s the overestimate 
is of order 10 km /s at [Fe/H] = —0.3). We observ e no depen¬ 
dence on [a/Fe[. Sanchez-Blazquez et al. (2011) emphasize 
that the effect is larger for young stellar populations, while 
our study here focusses on old stellar populations. Also, 
we are not interested in the velocity dispersions here, but 
in the parameters age, [Fe/H] and [a/Fe[. Any influence of 
this metallicity - velocity dispersion degeneracy would show 
up in our study as a bias in [Fe/H] recovery. Indeed, this 
may be the reason for the slight deviation from a perfect 
recovery at low metallicities observed in our simulations. 

We also explicitly test the influence of either normal¬ 
izing (rectifying) or not normalizing the spectra with a 
pseudo-continuum based on a running mean. We find that 
the fit with non-normalized continuum is formally slightly 
better, as expected due to the higher information content. 
We caution, however, that this is true for the ideal case we 
can test here, i.e. where the model and the "data" ideally 
match each other. Even a slight problem with flux cali¬ 
bration would potentially bias our results, we thus always 
normalize spectra from real data below. 

The age-metallicity degeneracy is a long-standing prob¬ 
lem in the study of old stellar populations. In the lower 
left panel in Figure [5 one can read of that the effect is 
of order O.ldex in [F^H] per 2 Gyr. Here we refrain from 
performing a formal fit to the simulated data. This would 
only pretend to provide more precise information as it is en¬ 
tirely unclear how the mock galaxy sample would relate to 
either the sample under study here or the rea l full g alaxy 
population. Comparison with Thomas et al. (2005 their 


that the trends are qualitatively the same, confirming once 
more that this effect is physical and does not depend on the 
fitting algorithm being applied. In particular it is surpris¬ 
ingly independent of whether the fitted stellar population 
is a SSP or the fitting algorithm allows for a more complex 
star formation history. To be more precise, the degeneracy 
is between [Fe/H] and age. There are degeneracies between 
age and [a/Fe], and between [Fe/H] and [a/Fe], but these 
are much less pronounced. 

A new problem in the analysis of integrated spectra from 
ETGs are the many lines of evide nce for a variation in the 
initial mass function of galaxi es ( van Dokkum fc Conroy 
2010b Cappellari et al. 2012). FerrAMateu et al. ( [2014^ 
explore how such variations would affect the analysis of 
the SFHs of ETGs and find that complete neglect of IMF 
variations would potentially significantly affect the results 
in particular in terms of the relative contribution of old and 
intermediate-age stars. Incorporating this caveat into our 
analysis is at present beyond the scope of this contribution. 

As described in Appendix [A] paradise performs a boot¬ 
strap estimate of errors on the recovered stellar population 
parameters. Figure|^verifies that these errorbars are appro¬ 
priate. Indeed, for correctly estimated Gaussian errorbars, 
the histogram of the quantity A(X)/ax should be a Gaus¬ 
sian with width 1, where A (A) = Ameasured ~ Atrue • This is 
almost exactly the case for age and the widths of the distri¬ 
butions in [Fe/H] and [a/Fe[. There are systematic offsets 
for [Fe/H] and [a/Fe[. Given the average error bars from the 
simulations these translate to a systematic overestimation 
of [Fe/H] by 0.07 dex and an underestimation of [a/Fe] by 
0.03 dex. Figure ^shows that the largest offset of --^0.1 dex 
occurs at low values of [Fe/H[. 


4.3. Can we resolve the enrichment histories? 

The term SEH implies the wish to represent the relation 
between age and abundance for each galaxy. To our knowl- 
edge this has not b een done in the literature. While [dej 
La Rosa et al.|(2011 ) have applied the spectral fitting tech- 
nique to the SFH, they have not actually resolved the SEH, 
preferring to quote a single value for the enhancement in 
a-elements. We thus here investigate in detail, whether our 
method is able to produce this kind of information reli¬ 
ably. We thus repeat the s ame plots of input vs. recov¬ 
ered quantities as in Section 4.2 for the properties [Fe/HL, 
[a/Fe[o, [Fe/H[i, [a/Fe],. The result is shown in Figure]^ 
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Fig. 6. Simulations of parameter recovery when splitting the 
star formation and enrichment history at a fiducial age of 9.5 
Gyrs. The dashed line in the lower right panel is a formal fit to 
the points in the data with a slope of ~ 0.3. 


To economize some space, we do not plot ages and metal- 
licities for the old and intermediate distinction. The reason 
is that the ages are fixed already by the way we separate 
the star formation history into bins and its recovered dis¬ 
tribution therefore carries no real meaning (although it can 
vary slightly within the width of that bin). Metallicity on 
the other hand is only a trivial combination of [Fe/H] and 
[of/Fe]. 

Figure |6 looks encouraging concerning the old stars in 
that clear relations scattering around the one-to-one line ex¬ 
ist, with an RMS width of about 0.1 dex. While it may not 
seem intuitive at first sight that the properties of the older 
stellar population, which contributes on average less to the 
total light than the younger or intermediate-age population, 
would be well-determined from a fit to the integrated spec- 
trum, this is actually in line with the findings of |Serra fc] 


Trager (2007). These authors found that the SSP-equivalent 
chemical composition depends mainly on the chemical com¬ 
position of the older of two constituting stellar populations. 
It is also in line with theoretical expectations that the ef- 


fects of [a/Fe 
stars ([Fe/H] 

(in particular Mg) are stronger for cooler 

Cassisi et al. 

2004 Coelho et al. 

2007 

San- 

som et al.||2013 1. 


This would in turn predict that the properties of the 
younger population are less well traced. Indeed, in the case 
of [Fe/H]i and [a/Fe]i the slope of the relation is clearly 
lower than one. A formal fit (the dashed line) yields slopes 
of ^ 0.3. To investigate the origin of this slope we have 
also investigated degeneracies between different quantities 
after this split into old and intermediate-age population. 
The same plots as in Figure for the separation of the 
SEH into old and intermediate populations shows a much 
larger scatter, but no noticeable increase in the degenera¬ 
cies, thus these plots are not shown here in the interest of 
space. However, and as shown in Figure the dominating 
degeneracies are those between [Fe/H], and [Fe/H]o and 
between [a/Fe]i and [a/Fejo- 

At this point we remind the reader of the classifica¬ 
tion into SLF and ACC galaxies introduced in Section 
Indeed, at specific points in this paper we separate galax¬ 



Fig. 7. The dominating degeneracies in the mean light-weighted 
values when separating SEHs into intermediate and old stars are 
those between the same properties of the two stellar generations. 
Each A is defined as the recovered minus the true value. 


ies in two types: I) Those where [Fe/H]o < [Fe/H];. For 
these the [Fe/H] evolution with time is consistent with self¬ 
enrichment by SNe la and these are therefore termed SLF- 
ETGs (represented by yellow dots in Figs. 17 to 22 below). 
2) Those galaxies where [Fe/H[o > [Fe/H];. Those are not 
consistent with self-enrichment and we infer that some ac¬ 
cretion of material must have taken place - whether these 
were mergers, gas accretion or anything else we remain ag¬ 
nostic about (such as winds red istributing elements within 
a galaxy, see Pipino et al.||2006 |. Nevertheless, for concise¬ 
ness we term those merger-accretion or ACC-ETGs (pink 
color in same figures). A more in-depth discussion of the 
physical usefulness of this distinction is given in Section 


We now explore how robust such a classification would 
be to degeneracies. Our input mock galaxy catalogue con¬ 
tains 50% SLF and 50% ACC galaxies. After fitting, 81% 
of the SLF galaxies would be correctly identified but only 
33% of the ACC galaxies. Furthermore, after fitting we 
would identify 74% of the sample as SLF galaxies and 26% 
as ACC galaxies. Finally, the number that is probably of 
most interest is the percentage of galaxies in each half of 
the parameters space that actually is what it pretends to 
be: The percentage of galaxies that would be classified as 
SLF galaxies after fitting and that really are SLF galaxies 
is 55%, while the percentage of galaxies that would be clas¬ 
sified as ACC after fitting and that are really ACC galaxies 
is 65%. While these numbers are sobering, one should re¬ 
member that they apply to a sample of mock galaxies that 
completely fills the parameter space available to our models 
and at a specified S/N of 40. This does not represent the 
real distribution of galaxies in this parameter space as we 
shall see below. In particular we shall explore in an obser- 
vationally motivated way how our conclusions depend on 
the exact location of the separation [Fe/H[o < [Fe/H[; + 
A, i.e. how this classification behaves within a range of 
-0.4 < A < 0.4. 

Besides direct classification we use another, more sta¬ 
tistically minded method to gauge the distribution of ob¬ 
served galaxies in the planes of [Fe/H[; vs. [Fe/H[o and 
[a/Fe[; vs. [a/Fe[o, namely the posterior distribution of fit¬ 
ted galaxies. We show the general case in Figure [^ one 
starting with a uniform distribution of mock galaxies in 
parameter space. For reasons to be seen later, another in¬ 
teresting case is the one where all galaxies start out in a 
narrow range of [Fe/H[; and [Fe/H[o abundances. We have 
chosen —0.2 < [Fe/H];, [Fe/H[o < 0.0. Clearly both cases 
produce significantly different posterior distributions in pa¬ 
rameter space that can be used to constrain the real dis- 
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Fig. 8. Comparison of the original distribution of parameters (upper row) vs. the one recovered after the fitting process (lower 
row). We show the two parameter spaces that are of most interest in the remainder of the paper, namely [Fe/H]i vs. [Fe/H]o and 
[a/Feji vs. [a/Fe]o. Left panel: The case in which the total library of mock galaxies is analyzed. Right panel: The case in which 
only mock galaxies with —0.2{ [Fe/H]i, [Fe/H]o (0.0 are analyzed. 


tribution of parameters to some extent - or at least can be 
used to exclude some real parameter distributions. 

In this section we have assessed the uncertainties intrin¬ 
sic to our method under the assumption that the model is 
a perfect representation of the data. We have found that 
these uncertainties exist, but that we can understand, de¬ 
scribe and quantify them. 


5. Fitting the observed spectra 

We now proceed to fit the sample described in Section 
with the method described in Section]^ We did not change 
the setup of the fitting in any way. This means that we 
fit the spectra with the same models, the same continuum 
rectification, and all other parameters being equal. There 
is only a slight difference in the spectra that may affect the 
fitting procedure to some extent. 

This is seen in the stacked fitting residual vectors of all 
galaxies as shown in Figure]^ Our sample has been selected 
to contain emission -line free galaxies as judged by |Brinch-| 
mann et al. (2004||. However, we find that the residuals 
do contain emission lines of H/3 and [OIII]. We attribute 
this finding to our excellent stellar population modeling 
and fitting procedure as well as to conservative estimates of 
emission line absence. In other words, while to the trained 
eye emission line residuals in stacked spectra may seem ob¬ 
vious, for an automatic detection procedure it is difficult 
to attribute a statistical significance to these features per 
galaxy, given similar-sized residuals in nearby regions of the 
spectrum. Emission lines affect only small regions of the 
spectrum and thus should not affect our overall results. In 
many cases they are masked away by our automatic outlier 
rejection in paradise (see Appendix [A|. We have verified 
that all our simulation results are inde]3endent of the po¬ 
tential missing information in the emission line region by 
applying an appropriate mask in wavelength space to the 
mock galaxies and rerunning all simulations. We found no 
difference in the results. 

As a side-note we mention the residual absorption fea¬ 
ture around 4700 A. Its origin is unclear, but it does not 
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affect our measurements as it is outside our fitted wave¬ 
length range. 

We also investigated the presence of any problems with 
data quality. Given the precision we are claiming to reach 
here (not much more than one percent in residual) we might 
be prone to data quality problems at a level not often seen 
otherwise. Problems might be present in the SDSS data, 
but also in the stellar templates. We note, however, that 
Figurej^does not indicate such problems. Indeed, the lower 
panel shows residuals in the observed frame where data 
problems such as sky emission line residuals should appear. 
The only obvious sky residual is at the very red end of 
the wavelength range shown, but way outside our fitted 
wavelength range. A paper dealing in much more detail 
with template mismatch vs. data quality is in preparation. 

Another way to assess our fit quality is to compare the 
cumulative histogram of obtained from the spectra in 
our data with the one expected from statistics, i.e. the 
distribution. This is done in Figure assuming 257 de¬ 
grees of freedom for our fitting procedure, as appropriate 
after taking into account the number of fitted pixels and the 
number of free parameters. The distributions do not match 
perfectly, but they are clearly very similar, both in shape 
and location. We believe this indicates that the fitting of 
SDSS spectra is indeed in the realm governed by x^ statis¬ 
tics, once the models are sophisticated enough and when 
restricting to a well-known wavelength range. 


In Sectionj^we used a velocity dispersion of 150 km/s for 
all tests, although the velocity dispersion was left a free pa¬ 
rameter in the fitting. Velocity dispersion is clearly different 
for every observed galaxy. We can, however, feed the SDSS 
pipeline velocity dispersion as a starting value to paradise, 
thereby essentially mimicking the fact that we know the ve¬ 
locity dispersion for the simulated galaxies. Indeed, we find 
that we do recover t he ve locity dispersion from the SDSS 
very well, see Figure liP 


^ There is one object out of 2286 where the SDSS pipeline ve¬ 
locity dispersion is 224.8 km/s, whereas paradise determines 
118.9 km/s. Visually this object clearly has been correctly fitted 
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Fig. 11. Comparison between the velocity dispersion from SDSS 
DR7 and those obtained in this work. 


Fig. 9. Mean residual vector and its standard deviation from 
all galaxies in the sample. The upper plot indicates the residuals 
in the rest frame of the galaxies, the yellow region has actually 
not been used to determine the best fit. [OIII] and H/? emission 
lines are labelled and clearly visible. The Mgb feature is also 
indicated for reference, but does not show evidence of specific 
mismatch. The lower plot shows the observed frame residuals 
indicating no major concern with the data. 



Fig. 10. Cumulative histogram of values obtained from the 
data (solid line) with the distribution expected from statistics 
(dotted line) for 257 degrees of freedom. 


udice to later or earlier results. They quote typical errors of 
1.48 Gyr in age, 0.04 dex in total metallicity, and 0.02 dex 
in [Fe/H] ratio, although it needs to be said for clarity that 
these earlier errorbars only take into account the formal 
uncertainties on the data. Our errorbars also include de- 
generacie s in the stellar populati on mix. In the more recent 
paper by Johansson et al. (2012) the uncertainty in [0/Fe] 
at S/N ~ 40 in the r-band is given as approximately 0.05 
dex, but it needs to be kept in mind that this is the uncer¬ 
tainty on a single element, while we group all a-elements 
together. 

Another important point to mention is that none of the 
derived galaxy ages are older than 13 Gyr, as we show in 
Section While this is predetermined from the range of 
model spectra in use, there is also no hint of an edge effect 
in the sense the galaxies that would be older than 13 Gyr 
cannot be fitted and therefore for those galaxies a higher 
metallicity would be fitted to compensate. We explicitly 
tested that this edge effect is present if we restrict our tem¬ 
plate library to SSPs with ages younger than 13Gyr. 


6. Results 

We now discuss our findings in terms of correlations be¬ 
tween the parameters we determined ourselves, i.e. ages, 
[Fe/H] and [a/Fe] and those available for our galaxies 
i.e. stellar mass M* and velocity dispersion a. 


The output of the fitting are tables with all parameters, 
such as mean light-weighted age ([Fe/H], [a/Fe]) of the en¬ 
tire spectrum, the same parameters only for the "old" stars 
(i.e. older than 9.5 Gyrs) and for the "intermediate-age" 
stars (i.e. younger than 9.5 Gyrs). Typical errorbars on our 
parameters are 0.2 Gyr on age and 0.01 dex on [a/Fe] and 
[Fe/H[. These errorbars reflect random errors and degenera¬ 
cies, but not model-related errors and systematic problems 
with the recovery. This can be compared with the typical 
errors quoted in one of the classical studies o f the field. For 
definiteness we choose Thomas et al. (2005), without prej- 


by paradise. We only have this one example, which however 
tentatively indicates that velocity dispersion mismatch would 
not have been a problem for our analysis. 


6.1. The a-a, M*-a and age-a relations 


In Figure 12 we plot [Fe/H] and [a/Fe] versus velocity dis¬ 
persion cr, stellar mass M* and mean light-weighted age. We 
quantify the presence or absence of correlations by means 
of Spearman rank coefficient test^ and by means of for¬ 
mal fits to each set of two parameters combinations using 


^ The Spearman coefficients measure the strength of a poten¬ 
tial correlation, where 0 is no correlation, 1 is a strong positive 
correlation and -1 is a strong anti-correlation. The probabilities 
on the other hand measure the signihcance of the correlation, 
i.e. the probability that the real correlation is zero. Therefore a 
value of 1 for this probability indicates that we know nothing 
about the real correlation, while a value of below 0.05 would 
imply that a correlation indeed exists at the the 95% confidence 
level. 
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Table 1. Correlation strengths for light-weighted mean param¬ 
eters 


Parameters 

Intercept 

Slope 

log(a) vs [Fe/H] 

0 . 02 ± 0.01 

-0.05±0.005 

log((T) VS [a/Fe] 

-0.46±0.01 

0.26±0.005 

log(M*/M 0 ) vs [Fe/H] 

-0.44±0.02 

0.03±0.002 

log(M*/M 0 ) vs [a/Fe] 

- 0 . 20 ± 0.02 

0.03±0.002 

age vs [Fe/H] 

0.15±0.003 

- 0 . 02 ± 0.001 

age vs [a/Fe] 

-0.09±0.002 

0 . 02 ± 0.001 



[Fe/H] 


Fig. 13. Correlation between [Fe/H] vs. [a/Fe] for the luminos¬ 
ity weighted average properties of the ETCs in our sample. The 
thick solid line is a line of constant, solar metallicity (Z=0.017). 
The two dashed line illustrate the extremes of the ETC locus, 
i.e. for [Z]=-0.15 and [Z]=0.18. 


6.2. Separating the old and intermediate components 

We now show results for resolved SEHs, i.e. looking at the 
abundance ratios of old and intermediate stellar popula¬ 
tions separately. Before we do so, as a cautionary remark, 
it is worth emphasizing that the uncertainty bars include 
the uncertainties introduced by the degeneracies character¬ 
izing galaxy evolution. Thus a direct visual comparison of 
the size of the error bar and the scatter in the figures is 
misleading. The size of the error bar indeed denotes the 
range of possible values for each single galaxy. However, on 
average, galaxies will tend to be assigned the mean value 
of that possible range. Scatter in figures is mainly driven 
by real observational uncertainties and does not represent 
the degeneracy uncertainties as assessed by paradise. That 
said, mean uncertainties on our parameters are 0.17 dex for 
[Fe/H]i, 0.14 dex for [Fe/H]o, 0.12 for [a/Fe]i, and 0.13 for 
[a/Fe]o. 



the LINFIT module in IDL. The results of the fits are given 
in Table We confirm correlations between log(cr) as well 
as log( MVM ( 7 )) and [a/Fe] as found in previous work (see 
Section [O] ) . In Figure [I^ we also see a strong correlation 
between age and [a/Fe] as well as [Fe/H] in the sense that 
[a/Fe] decreases for younger light weighted mean age of the 
galaxy, while [Fe/H] increases. 


There seem to be onl y weak correlations be tween 
[Fe/H] and log/cr). However, Johansson et al.| ( 2012| | show 
that when the sample is decomposed into age bins, each age 
bin by itself shows a strong correlation with log(cr). In the 
interest of space we refrain from repeating the exercise. 


Additionally the relation between age and [a/Fe] shows 
a clear change of slope at age ^9 Gyr. We have fitted the 
two regimes separately and find that [a/Fe] — -0.01±0.004 
+ O.OliO.001-age for galaxies younger than 9 Gyr, while 
[a/Fe] - -0.20±0.007 + 0.03±0.001-age for galaxies older 
than 9 Gyr. For both regimes the probability of the absence 
of any correlation is zero. 

The two rightmost panels of Figure [^already indicate a 
correlation between [a/Fe] and [Fe/H], which we explicitly 
plot in Figure [T^ We find that early-type galaxies all have 
on average a total metallicity that is roughly solar. It is 
the ratio between the abundances of Fe and the a-elements 
that changes between galaxies. 
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Fig. 14. Average star formation history of all galaxies in the ob¬ 
served sample in terms of the present day contribution of each 
stellar population to the total luminosity of the galaxy. All galax¬ 
ies contribute equally to this average, i.e. the total luminosity 
for each galaxy has been normalized to one before averaging. 


First we repeat in Figure p^ the upper panel of Figure 
[2 i.e. we show the non-binned average star formation his¬ 
tory of all galaxies in our sample. As this is the observable, 
we show the star formation history in terms of the lumi¬ 
nosity contributions of stars of each specified age to the 
total luminosity of each galaxy. It is noteworthy that sev¬ 
eral peaks in the star formation history appear, namely at 
13 Gyr, 9 Gyr and a somewhat broader one around 6 Gyr. 
Our SSPs are sampled in steps of IGyr, therefore interme¬ 
diate ages could in principle be assigned significant weights 
by paradise. On one hand it is entirely possible that this is 
an artifact of our fitting procedure, given that such features 
do not appear in the overall star formation history of the 
universe (Lilly-Madau plot). On the other hand, the two 
older episodes of star formation seem entirely compatible 
with the classical two- infall scenario of chemical evolution 
(Ghiappini et al. ]T9^ that was originally proposed for the 
Milky Way. Note that these two oldest peaks also contribute 
to justifying our separation into old and intermediate-age 
stars at a lookback time of 9.5 Gyr. We come back to this 
point below when discussing SLF vs. AGG-FTG percent¬ 
ages. 
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Fig. 12. Correlations of [Fe/H] and [a/Fe] with velocity dispersion, stellar mass and mean light-weighted age for our sample 
galaxies. It is clear that the tightest relations are those with age. Note also the change of slope in the relation between age and 
[Q:/Fe]~9Gyr. Average error bars are given in the upper right corner of each panel and spearman rank correlation coefficients are 
also noted. 


Table 2. Correlation strengths old stars only 


Parameters 

Intercept 

Slope 

log(cr) vs [Fe/H[o 

-0.73±0.046 

0.29±0.020 

log(cr) vs [a/Fe[o 

0.37±0.046 

-0.05±0.020 

log(MVM 0 ) vs [Fe/H[o 

-0.83±0.066 

0.07±0.006 

log(M*/M 0 ) vs [a/Fe[o 

0.46±0.063 

-0.02±0.006 

age vs [Fe/H[o 

0.13±0.013 

- 0 . 02 ± 0.001 

age vs [a/Fe[o 

0.39±0.012 

-O.OliO.OOl 


Finally, there is a small peak at less than 1 Gyr which is 
outside the coverage of our [a/Fe]-enhanced SSP templates. 
As described in Section we include five templates with 
younger ages, but with solar element ratios only. The total 
contribution to the light of the population younger than 2 
Gyr is however less than 5%. Due to its lower mass to light 
ratio than the older stellar populations this translates to an 
even smaller contribution to the total mass of the galaxy. 
We thus confirm that ETGs harbor very little young stars 
maybe even none, if currently unmodeled old and hot 


stellar populations indeed do play a role (see e.g. Ocvirk 


20101 . 


It is interesting to start the analysis with plotting the 
properties of the old stellar populations vs. the global pa¬ 
rameters of their host galaxies. We refrain from doing so for 
the intermediate-age populations due to the strong degen¬ 
eracies we identified in Section |4.3| and the lower recovery 
success for the intermediate component with respect to the 
old one we identified there. Also, our sample is dominated 


by galaxies in which the old stellar population contributes 
more than half of the total light, contributing to possible 
uncertainties in the properties of the intermediate-age stel¬ 
lar populations. Figure [T5| shows that the correlations seen 
for the global stellar populations do not hold up for the 
older stars in the galaxies (compare also Table 1^. Indeed, 
while the Spearman rank test indicates a very low prob¬ 
ability for the dependence of [a/Fe]o and [Fe/H]o on age 
to be spurious, inspection of the plot tells another story. 
There is a cloud of points at older ages, while the galaxies 
with younger ages show much larger scatter. This can be 
understood as being caused by difficulties in accurately es¬ 
timating [a/Fe]o and [Fe/H]o for galaxies with significant 
intermediate-age populations. We therefore attach no im¬ 
portance to these relations. For a and M* vs. [a/Fe]o, the 
Spearman test indicates weak correlation and indicates a 
high probability (0.68) of no correlation for M* vs. [a/Fe]o. 
We are thus left with the indication that [Fe/H]o depends 
positively on mass as parameterized by either tr or M*, 
but more strongly on the former. This correlation is driven 
mainly by the presence of objects with low [Fe/H]o at low 
a values. Overall these trends are opposite to the aver¬ 
age properties, where [a/Fe] depends clearly on mass and 
[Fe/H] does not. 

The two rightmost panels of Figure [^already indicate 
that the correlation between [a/Fe] and [Fh/H] that is valid 
on average for early-type galaxies does not hold up when 
considering the old stars only. Indeed, in Figure [16] we find 
that old stars in early-type galaxies have on average super¬ 
solar metallicity, where both [Fe/H[o and [a/Fe[o are super- 
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Fig. 15. Correlations of the [Fe/H]o and [a?/Fe]o with the global properties of the sample galaxies a, M*, and age. Average error 
bars are given in the upper right corner of each panel and spearman rank correlation coefficients are also noted. 



old 


[Fe/H] 


Fig. 16. Relation between [Fe/H]o vs. [a/Fejo, the thick line is 
a line of constant, solar metallicity (Z=0.017). This is the same 
as Figure but only for the stars older than 9.5 Gyr in each 
galaxy. The old stars in early-type galaxies, have on average high 
total metallicity, significantly super-solar. Average error bars are 
shown in the upper right corner. 


solar. That these properties seem to be independent of mass 
(see Figure 151 might be one of the most puzzling findings 
of this paper. 


6.3. Exploring the possible dichotomy SLF vs. ACC 

In Figure we show the distribution of [Fe/H] and 
[a/Fe] of the intermediate vs. the old component (left pan¬ 
els). In [Fe/H] galaxies can lie both above the identity line 
and below it. In [a/Fe], almost all galaxies lie below the 
identity line, with a few slight outliers that can be thought 
of as being due to uncertainties. When looking back at Fig¬ 
ure]^ and the discussion in Section [T3l it is obvious that any 
statement on the distribution in these panels must take into 
account the uncertainties due to the degeneracies in prop¬ 
erties between old and intermediate-age stars. In turns out 
it is easier to discuss the results in terms of what the data 
are incompatible with and this is what we do. 

The distribution of points in the upper left panel on 
Figure is incompatible with a flat distribution in that 
plane. I.e. early-type galaxies do not have a random dis¬ 
tribution of [Fe/H[i and [Fe/H[o. Rather their distribution 
is more similar to the right panels in Figure [^ where the 
underlying sample had a narrow distribution of[Fe/H[i and 
[Fe/H[o centered around -0.1 for both. 

The distribution of points in the lower left panel of Fig¬ 
ure [^ is equally incompatible with a flat distribution in 
that plane. I.e. early-type galaxies do not show a random 
distribution of [a/Fe[i and [a/Fe[o. It seems rather secure 
that the oldest generation of stars has a high [a/Fe] ratio, 
higher than 0.1, whereas the intermediate-age stars have 
a lower ratio of [a/Fe], certainly below 0.1 and going as 
low as our model permits, i.e. -0.2. Note that [a/Fe] values 
below zero for late stellar generations are a generic pre¬ 
diction of chemical evolution models in massive galaxies 
([Fe/H] [Minchev et al. 


2013 their Fig. 5). On the other 


hand low-mass galaxies do not necessarily show low [a/Fe[. 
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Fig. 17. Left panels: The distribution of mean light-weighted [Fe/H] and [a/Fe] for the old vs. intermediate sub-populations. 
Galaxies in yellow are those where [Fe/H]o < [Fe/H]i-0.15 (SLF-ETGs), pink where [Fe/H]o > [Fe/H]i+0.15 (ACG-ETGs). Galaxies 
with uncertain classification are shown in black. Average uncertainties including all degeneracies are shown in the upper left corners. 
Right panels: The age dependence of [Fe/H] and [a/Fe] within each galaxy, i.e. their enrichment history. Galaxies with uncertain 
classification have been omitted for clarity. 


Recchi et al. (20011 for example show that stellar winds 
tend to carry elements produced by SNela away, thus leav¬ 
ing behind gas with higher [a/Fe] than expected from the 
extent of their SFH. 


Coming back to the distinction discussed in Section [T^ 
between ETGs that are compatible with self-enrichment 
(SLF-ETGs) and mergers/accretion (ACG-ETGs), we wish 
to identify galaxies that we can reasonable safely assign 
to any of these categories. SLF-ETGs are those where 
[Fe/H[o < [Fe/H[i - A, and ACC-ETGs are those where 
[Fe/H[o > [Fe/H[i + A. Motivated by the sizes of the er¬ 
ror bars quoted in the beginning of Section |6.2| and after 
some testing in the simulations, for definiteness we choose 
A = 0.15. Our results qualitatively do not depend on the 
exact size of A. We emphasize that this choice of A also 


For this A and for the simulations, the percentage of galaxies 
identihed as ACC-ETGs after the fitting process, which are truly 
ACC-ETGs is 77%. The percentage of galaxies identihed as SLF- 
ETGs after the Htting process, which are truly SLF-ETGs is 
63%. 


mitigates concerns over the su btle drop in metallicit y at 
late times shown in the paper by Vazdekis et al. (19961 and 
discussed in the introduction. We call the remaining galax¬ 
ies that are neither one nor the other the ’gray-ETGs’. Note 
that for the observational results we include into the class 
of ’gray’ galaxies those galaxies that have less than 15% of 
light in either the old or intermediate stellar populations. 
This is to account for possible uncertainties in stellar pop¬ 
ulation modeling and in the fitting algorithm. 


For the sample of 2286 galaxies analyzed here we find 
that 29% belong to the AGG-ETG class. 63% are unclassifi- 
able, but most of those would be called AGG-ETGs as well if 
not for the requirement to have 15% of light at least in the 
intermediate-age population. Only 8% of the galaxies are 
bona-fide SLF-ETGs, i.e. consistent with a self-enriching 
star formation history. 


We caution that effects of degeneracies between the 
properties of the intermediate and old stellar popula¬ 
tions, sample selection effects, the uncertainty on the exact 
boundary to be applied between old and intermediate stel- 


Article number, page 13 of 


22 


















A&zA proofs: manuscript no. sdssetgs 


lar populations, specific effects of chemical evolution such 
as the drop in metallicity for young stars in massive ETGs, 
all contribute to make these numbers indicative only. To 
indicate the possible size of changes due to at least one 
of these effects we have recomputed the percentages of SLF 
and ACC-ETGs for a choice of 8.5 Gyr as the boundary be¬ 
tween old and intermediate-age ETGs. We find that this in¬ 
verts their relative importance, i.e. 30% of the sample would 
then be SLF ETGs, while only 3% bona-fide AGG-ETGs re¬ 
main. This high sensitivity of the results to the adopted age 
boundary is due to a very low [Fe/H] abundance attributed 
by the fitting code to the spike in the SFH at 9 Gyr. This 
significantly lowers the average [Fe/H]o abundance, thus in¬ 
creasing the number of SLF-ETGs. 

Given the magnitude of this effect, we have to care¬ 
fully state our arguments for choosing 9.5 Gyr lookback 
time as our boundary. Indeed, we have argued in the intro¬ 
duction on physical grounds that 9.5 Gyr is a more physi¬ 
cal boundary. On technical grounds, we refer to Figure 
which shows that the 9 Gyr age bin in the finely sampled 
SFH has a close significant neighboring age bin at younger 
age (7 Gyr), while the next significant old age bin (at 13 
Gyr) is more distant in look-back time. Generally speaking, 
it is clear that there are significant degeneracies between 
the properties of adjacent bins in the reconstructed SFH. 
We have shown the extent of these degeneracies in Figure 
[^for the coarsest possible choice of SFH, namely an SFH 
with two age bins. Degeneracies only increase for an SEH 
with a finer age grid. Degeneracies between the 7 and 9 Gyr 
age bins shown in Figure are therefore much more likely 
than between the 9 and 13 Gyr age bins. This is the tech¬ 
nical reason why we believe the results based on a 9.5 Gyr 
boundary to be more reliable and significant than those we 
would have obtained based on any other boundary. 

Figure [T8| repeats Figure but adding the information 
about the enrichment history classification. It is immedi¬ 
ately clear that on average and for our sample AGG-ETGs 
are significantly older than SLF-ETGs (compare also Figure 
25). Also, on average SLF-ETGs have lower velocity disper¬ 
sion and lower [a/Fe] than AGG-ETGs. We come back to 
these features in the discussion Section [H 


6.4. Environment 

While nearly out of the focus of this paper, environment 
is such an important parameter for galaxy evolution that 
it would be foolish to leave it entirely out. Much work has 
been done on assessing the correlations between environ¬ 
ment and the properties of early-type galaxies and we do 
not wish to simply repeat this work here, as we strongly 
expect to confirm it. Rather, in the interest of space, we 
concentrate only on the correlation of the old stellar pop¬ 
ulation properties [Fe/H]o and [a/Fe]o with environment, 
taking into account our classification in AGG-ETGs and 
SLF-ETGs. We are particularly interested to see whether 
the initial star formation eve nt indeed depends on environ¬ 
ment ( [Fe/H] Renzini|[2006 1. 

As a measure of env ironmental density we use the one 
by Tempel et al. (|2012 1, specifically we use the density as 
derived within 8 Mpc, i.e. their parameter dens8. Their den¬ 
sities are based on luminosity densities around each galaxy 
calculated in three coordinates, i.e. two projected spatial 
distances and one based on redshift, the latter however in¬ 
cluding advanced techniques to suppress peculiar motions 
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Fig. 19. Environmental density vs. light-weighted [Fe/H] and 
[a/Fe] (upper panels) and vs. [Fe/H]o and [a/Fe]o of only the 
old stars in each galaxy. The mean for each subpopulation is 
also shown as a large square with black borders. 


such as suppression of the finger-of-god effect. We refer to 
the original paper for all the details. 

Perhaps surprisingly we find in Figure [T^ that the old 
stellar population may not be affected by environment. This 
is clearly differ ent from what other authors fo und for the 
entire galaxies ( Sanchez-Blazque z et al.|| 2006a|). It h owever 
agrees with inferences made in [Thomas et al.| ( |2010 l). These 
authors discuss that ETG formation has occurred in two 
phases, where ETG galaxy f ormation depended o n mass 
only and not on environment. Rettura et al. (20101 on the 
other hand argue that while the environment may modulate 
the timescale of the star formation history for ETGs, the 
formation epoch is essentially independent of environment 
and seems to depend on mass only. We here add another 
piece to the puzzle in saying that for the oldest stars in each 
galaxy, older than 9.5 Gyr, galaxy formation seems to have 
been independent of mass and environment, at least for the 
mass range probed here. 


7. Aperture effects 

As discussed in Section the cores and the envelopes of 
ETGs may have assembled in different ways and at dif¬ 
ferent times. Alsi^winds can redistribute elements within 


a galaxy (Pipino et al. 20091. Observationally, gradients 


([Fe/I 
for Ef 

1] |Sanchez et al. 

2013 

Perez et al. 

2013 

|. However, 

L'Gs the gradients in age and [a/h'e 

are observed to 

be small or null on average 

(IMehlert et al. 2003 Rawle 

et al. 

2008 

Kuntschner et ai.| 20101. Not so for the era- 


dient in |Fe/H|, which can be significant ([Fe/H] 


et al. 2008|). The SDSS uses a single fiber to observe the 


centers of the target galaxies. The fibers fixed aperture on 
the sky (3") covers a different fraction of the target galaxy, 
depending on its intrinsic size and redshift. This effect is 
somewhat counteracted by the effects of seeing, as seeing 
will also scatter light from larger radii into the fiber area 
in the focal plane. It is nevertheless clearly important to 
assess this covering fraction and the potential impact on 
ou r results. Proper ape rture corrections as demonstrated 
Gallazzi et al. (2008) is beyond the scope of the present 


paper. We also refer to the discussion in Ghoi et al. 


(20141 
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Fig. 18. Correlations o f [Fe /H] and [a/Fe] with velocity dispersion, stella r m ass and mean light-weighted age for our sample 
galaxies (same as Figure 121. Colors reproduce the classification from Figure 17 into SLF-ETGs (yellow) and ACC-ETGs (pink). 
The mean for each subpopulation is also shown as a large square with black borders. Galaxies that could not be reliably classified 
are shown as black dots. 


for another angle on the same question, albeit they work 
on a larger redshift range. 

Figure [^shows the dependence of key physical proper¬ 
ties on redshift for our sample. In the uppermost panel the 
effects of apparent magnitude selection are obvious in that 
we sample more massive objects at larger distances from 
Apache Point Observatory. The next panel shows the de¬ 
pendence of the covering fraction on redshift in terms of the 
ratio between the fiber size and the galaxies effective radius. 
We use the effective radius from the SDSS pipeline surface 
brightness profile fit with a de Vaucouleurs profile. It illus¬ 
trates that two effects nearly cancel each other: on the one 
hand, the higher the redshift, the larger the area covered 
by a fiber is physically smaller on the galaxy. On the other 
hand, the galaxy itself will tend to be more massive and 
therefore larger. The covering fraction thus depends less 
on redshift than the stellar mass. The next panel shows the 
expected non-dependence of [a/Fe] on redshift, given that 
[a/Fe] varies little with mass and galaxies have flat gra¬ 
dients if any. The final panel shows that [Fe/H] depends 
in an interesting way on redshift. Because of the stronger 
gradients and stronger dependence of [Fe/H] on mass it is 
nontrivial and beyond the scope of the present paper to 
understand this. 

The crucial test, however, is shown in Figure [^ 
where we show the dependence of [Fe/H] on covering frac¬ 
tion at a fixed mass (here chosen to be 10.5<log(stellar 
mass/M 0 )<ll.O). There is no obvious bias to be seen, be¬ 
sides a potentially spurious tendency to have less [Fe/H[- 
rich objects at larger covering fraction. We conclude that 
tests on possible aperture biases are not conclusive enough 
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Fig. 20. The dependence of key physical properties of our 
sample galaxies on redshift. From top to bottom: stellar mass, 
[a/Fe], age and ratio between radius covered by fiber and ef¬ 
fective radi us o f the galaxy. Golors reproduce the classification 
from Figure [T t] into SLF-ETGs (yellow) and AGG-ETGs (pink). 


to endanger any of the fairly general conclusions of this 
paper. 


Any further investigation of aperture effects, in partic¬ 
ular as concerning the resolved SEHs, i.e. the potentially 
different biases on the old and intermediate-age stars) will 
have to await analysis of the currently on going integral 


field spectroscopy surveys such as CALIFA (Sanchez et al. 
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Fig. 21. Dependence of [Fe/H] on coverage fraction at fixed 
mass. No b ias is obvious. Colors reproduce the classification from 
Figureinto SLF-ETGs (yellow) and ACC-ETGs (pink). 


2012|, SAMI dAllen et al.||2015|), and Manga (IBundy et al. 

WEi -- - 


Fig. 22. Locus of our s ample galaxies as com pared to the mass- 

I05|l (solid line, dotted 
Delgado et al. 


metallicity relations of Gallazzi et al. 
lines are 1 a scatter) and of Gonza 


(20141 


(dashed line). While there is general agreement, the clear differ¬ 
ence in slope within the mass range of our galaxies is interesting 
and may point to sample selection effects. 


8. Discussion 

8.1. Comparison of average quantities to earlier observational 
work 

A simple test is to situate our galaxies on the stellar mass 
- stell ar m etallicity relation as defined b y [Gallazzi et al.| 
(20051 and Gonzalez Delgado et al. (20141 as done in Fig¬ 
ure]^ The overall match is satisfactory, given the entirely 
different analysis methods. Note though that we use the 
same data as for the Gallazzi relation (although our sample 
is smaller) and we use their stellar masses. Any difference 
between the two is therefore directly tied to the derivation 
of the metallicities and abundance ratios. Of particular in¬ 


terest is the d ifference in slope. Both the Gnllazzi et al.j 
(20051 and the Gonzalez Delgado et al. (20141 relation im- 
ply a significant slope at masses around 10 ^^ Mq , where 
our analysis seems to imply n one. This slope i s conf irmed 
for the case of ETGs alone in Gallazzi et al. (20061. The 


absence of a slope in our analysis may point to a sample 
selection effect at the low mass end of our sample or may 
be related to the overall offset in total metallicity. 

Another interesting difference to earlier studies is the 
range covered by the [a/Fe] parameter. None of our sam¬ 
ple galaxies scatters above 0.3 in its light-weighted average 

[g/Fe] 


earlier work ([Fe/H] Thomas et al.|2005 

Kuntschner et al. 

20101 Johansson et al. 2012|) may be due to observational 

scatter, for the stacked spectra of 

Graves et al. 

([ 20101 , this 


lutionary effects of [g/Fe] into account and simplified the 
corrections for the spectral effects. They would thus under¬ 
estimate the Mg feature at give n g-enhancement (Goelho 
et al. ]|2M7l jVazdekis et al.|[2011[ ). Really pinning down the 
cause of the offset is out of the scope of the present paper 
though. 

There are many works in the literature r eporting corre¬ 
lations between [g/Fe] and either cr or M*. Thomas et al. 
(2005| find slopes of 0.28 for [g/Fe] vs. log(CT) and 0.06 for 
I g/Fe] vs. log(M*/M 0 ). Trager et al. ( 2000| ) on the other 
hand find a steeper slope of 0.33 for [g/Fe] vs. log((T), simi¬ 


lar to Nelan et al. (|2005| who q uote 0.31 ± 0.06 for the same 
slope. Finally, Hernardi et al. (20061 quote a slope of 0.32, 
while Graves et al. (2007|) give 0.36. Our slope of 0.26 seems 
to lie at the low end of the scatter in the literature, which 
is probably tied to the lower maximum [g/Fe] reached by 
our analysis. 

A correlation that is more rarely plotted is the one of 
[g/Fe] and [Fe/H] with age (Fig. 12 right panels). Gor- 
relatio ns of age and [g/Fe] were aEiown in Graves et al. 
(2 0 lot from stacked s pectra (their Fig. 4)_and per galaxy 
Kuntschner et al. (20101 (their Fig. 6 )j^ In both cases 


the change of slope at an approximate age of 9 Gyr is visi¬ 
ble, while the overall normalization in [g/Fe] seems shifted 
higher as compared to our results by 0.05 to 0.1 dex. Due 
to the nature of the stacked spectra in [Graves et al. (|2010l 


and the small sample size in Kuntschner et al. 


we 


cannot comment on a comparison of the scatter around the 
mean relation. Although a direct comparison is difficult be¬ 
cause of the di fferent [g/Fe] definitio n, a change of slope 
was also seen in Gallazzi et al. (20061, albeit at ages of ^6 
Gyr, instead of OGyr.We also note that these re sults are in 
contrast to the earlier work of jjprgensen (19991 who found 
no correlation between [g/Fe] and age. 

For the same correlation, it i s important to heed the 
warnings of Thomas et al. (2005[|, who discuss the impor¬ 
tance of degeneracies when using age as a parameter. They 
find that for index analyses [Fe/H] and [g/Fe] are anti¬ 
correlated with age. As our dependency of [g/Fe] on age 
is positively correlated with age, the observed slope cannot 
be due to degeneracy. For [Fe/H] we find that the slope of 
the degeneracy is approximately -0.05 dex/Gyr (see Fig.|^, 
whereas the slope of the relation is -0.02 dex/Gyr. Further¬ 
more our relation extends over 10 Gyr, which is a much 
larger range in age than what would arise if all galaxies 
had the same age and would extend on an age-[Fe/H] re¬ 
lation purely due to a degeneracy. We thus conclude that 


de La Rosa et al. 


(20111 have shown a correlation between 


the percentage ol mass in old stars and [a/Fe], but that is more 
difficult to compare to our results. 
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Fig. 23. Comparison of the old and intermediate-age stellar 
populations in ETGs with the properties of the stars in the 
Milky Way. There is a fundamental difference in the fact that 
for ETGs we measure the luminosity weighted average of the 
properties over all stars in the galaxy. For the Milky Way on the 
other hand, the averages are over single stars and not all types 
of stars contribute to the average. 


the age-metallicity degeneracy does not give rise to the ob¬ 
served correlation. 

These correlations of [Fe/H] and [g/Fe] w i th age ex¬ 
plain the findings by Sanchez-Blazquez et al. (2006a|), in 
which metallicity as measured by different indices (Fe, 
MgB) seemed to behave differently as a function of age. Fe- 
sensitive indices are a reasonable proxy of [Fe/H], whereas 
the MgB indices depend on age, [a/Fe] and [Fe/H], thus 
likely washing out the pure [a/Fe] signal as a function of 
age. 

We have shown in this paper that the age - [a/Fe] re¬ 
lation can be p ut on a firm found ation using the advances 
presented here. Yates et al. (20131 show that even the stel¬ 
lar mass - [a/Fe] relation alone is a good constraint on the 
delay time distribution of SNe la. We suggest that the age 
- [a/Fe] relation should provide additional interesting con¬ 
straints, especially when considering the two regimes with 
different slopes. 

Fi nally, w e caution that results such as those in [Trager] 


et al. (20081 seem to indicate that while many of the cor¬ 
relations we discuss are true on average over large samples, 
every single structure or galaxy cluster may have its own 
peculiar history. 


8.2. Comparison between SLF-ETGs and the Milky Way 

We have been able to diagnose for the first time the physical 
properties of old and intermediate stellar populations of 
ETGs. This means that for the first time we can compare 
the properties of these sub-populations to measurements 
otherwise only available for very nearby galaxies, mostly 
within the local group. We attempt this in Figur e]^ which 
reproduces the cartoonish expectations of Figure [ijbut this 
time with real data. 

The two systems in the Milky Way that can be com¬ 
pared to ETGs in a more meaningful way are the thick 
disk above the plane (where contamination from the thin 
disk is minimal) and the bulge. For the thick disk w e plot 
the fiducial line [a/Fe] = —0.2x[Fe/H[ -1-0.10 from Nide- 



[Fe/H] 



[Fe/H] 


Fig. 24. Old and intermediate age stellar populations in ETGs 
separated by the classification according to their enrichment 
history. SLF-ETGs (upper panels) are consistent with self¬ 
enrichment, while ACG-ETGs must have suffered either strong 
merging/accretion or strong internal winds. The properties of 
the merged galaxies can be tentatively read of as the properties 
of the intermediate age population, i.e. low in [a/Fe] and lower 
in [Fe/H] than the main body of the galaxy. 


ver et al. (20141. We also add the Magnesium abundances 
of the bulge stars o lder and younger than 9.5 Gyr from 
Bensby et al. (20131. These represent bona-fide bulge stars 
and due to the reliable ages from micro-lensing we are able 
to apply the same selection as for the old stars in our ETGs. 
For clarity we have not plotted any values for the thin disk 
of the Milky Way, which shows continuous star formation 
until today, ft would lie between [Fe/H] = -1 and 0 and 
around 0.0 in [a/Fe[. 

When looking at Figure[^it is important to appreciate 
the difference in the probed stars. For the ETGs analyzed in 
this work each red dot represents the luminosity weighted 
average abundance values of all stars older than 9.5 Gyr in 
the galaxy. For the bulge and the thick disk, each star for 
which abundances have been measured contributes equally 
to the locus that is reported in this figure, i.e. we compute 
a simple mean. We proceed to formulate some basic con¬ 
clusions that are robust against possible selection effects, 
given the very different methods to obtain the data in MW 
and ETGs. 

ft is immediately clear from the figure that normal, mas¬ 
sive ETGs differ wildly in their chemical enrichment history 
from both the thick disk of the Milky Way and the bulge of 
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the Milky Way. Obviously the initial star formation event in 
ETGs must have been much more efficient in forming stars 
than even the stars of the thick disk or th e bulge. Our state- 
ment is complementary to the findings of Tang et al. (20141 
who find that the bulge is closer to ETGs than normal disk 
stars in its enrichment history. However, here we empha¬ 
size t hat the way in wh ich enrichment proceeds is different, 
while Tang et al. (20141 find that the result from the enrich¬ 
ment (re7T[ie"3Ifferent element abundances) is different. To 
gain a bit more insight into the variety of manifestations of 
this fact in our sample of ETGs, we reproduce Figure 
for the two classes of ETGs defined in the present paper, 
namely SLF-ETGs and AGG-ETGs. In Figure we are 
particularly interested in the SLF-ETGs which at least are 
compatible with self-enrichment and a causally connected 
stellar enrichment history. Indeed, we find for these galax¬ 
ies that the intermediate age stellar populations follow a 
relation between [Fe/H]; and [a/Fe]i that is reminiscent of 
those followed by the thick disk and bulge stars, albeit with 
a steeper slope and extending to lower values of [a/Fe]i. 
Very low values of [a/Fe] are indeed expected for the stars 
formed afte r a long enrichment history (see [Fe/H] Minchev 
et al ||2013 their Fig. 5). 

IT indeed, two different processes (self-enrichment 
vs. merger/accretion/winds) leave their imprint on the 
abundances of the intermediate age stars in ETGs, then 
it is to be expected that both processes act in all galax¬ 
ies. SLF-ETGs and AGG-ETGs would then mostly differ 
in which of these two processes dominate, while the other 
would still play a role. This might be the explanation for 
the steeper gradient in the [a/Fe]i vs. [Fe/H]i relation as 
compared to the bulge stars. 

While such comparisons between wildly different galax¬ 
ies that have been analyzed in wildly different ways should 
be taken with caution, our results are also consistent with 


those of Sil’chenko (1993); Proctor & Sansom (2002). Both 
of these author groups nnd that the bulges of disk galax¬ 
ies are on average younger and less enhanced in g-ele ments 
than ETGs. In particular, Proctor & Sansom (2002) state 
that spiral bulges show a positive correlation between ve¬ 
locity dispersion and metallicity, as expected from a sim¬ 
ple scenario in which mass builds up slowly over a Hubble 
time, whereas ETGs show an anti-correlation as also found 
in our data. That our data do not support the presence 
of a clearly identifiable knee as we are used to from Milky 
Way data is consistent with this finding of strongly diverg¬ 
ing processes in the chemical evolution history of late and 
early type galaxies. 

Finally, we not e that we are unable t o directly confirm 
the predictions by Lackner et al. (2012), as in these sim¬ 
ulations the accreted stars are ~2.5 Gyr older, and ~0.15 
dex more metal poor on average than the stars formed in 
situ. This does not seem to be compatible with the distri¬ 
bution in [F e/H| vs. [a/Fe| space we descri be here. On the 
other hand, Navarro-Gonzalez et al. (2013) point out that 
it is exactly these older and metal-poor stars that are in the 
outskirts of galaxies, thus providing a potential way out, as 
our SDSS spectra overweight the inner parts of galaxies. 


8.3. On prolonged star formation in SLF-ETGs and the 
properties of minor mergers in ACC-ETGs 


Figure 25 shows the star formation histories (i.e. the SEHs 
marginalized over all abundances) for the three different 
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Fig. 25. Average star formation history of galaxies in the three 
classes as classified by their enrichment history. The figure shows 
the present day contribution of each stellar population to the 
total luminosity of the galaxy. All galaxies contribute equally to 
the average within each class, i.e. the total luminosity for each 
galaxy has been normalized to one before averaging. 


classes of galaxies introduced in this paper. We remind the 
reader that there are some caveats to this classification in 
terms of the exact location of the age boundary used to 
separate old from intermediate age stars. The following dis¬ 
cussion is, however, insensitive to these caveats. A general 
feature of Figure is that the percentage of intermediate 
age stars is much higher in SLF-ETGs than in AGG-ETGs. 
This would seem to indicate that mass growth through late 
in-situ star formation can be significant. On the other hand, 
late mergers only very rarely seem to contribute sufficient 
mass to significantly change the overall SFH of the ETGs. 
This behavior has been called quasi- monolithic (or early 
hierarchical) scenario in the literature ( Merlin et al.||2012 ). 
For the moment this inference does not contradict a picture 
in which significant late mass growth occurs through minor 
mergers, as these are expected to be deposited at larger 
radii, outside the coverage of the SDSS fiber. Also, this 
finding is consistent with the recent results of |Hahn et al.| 
(2014) in which the authors find from entirely different lines 
of evidence that most of the ETGs stopped forming stars 
due to internal processes and not due to external influences. 
Our dat a, however, are i n ten sion with simulations such as 
those by Lackner et al. (2012), in which the authors found 
that minor mergers contribute between 15% and 40% of 
the total mass and that the percentage of mass brought in 
by minor mergers increases with stellar mass. We do not 
find that accretion increases with stellar mass, but rather 
that in-situ, prolonged star formation decreases with stellar 
mass. 


If, indeed, the intermediate age stars in AGG-ETGs 
mostly come from merging over a Hubble time, the proper¬ 
ties of these stars should carry information on the proper¬ 
ties of the galaxies in which they formed, i.e. on the prop¬ 
erties of the minor mergers before the merging. In Figure 
|24]the fit results for those stars that contribute to the in¬ 
termediate age population clearly tend towards the lower 
boundaries imposed by the model in terms of [Fe/H] and 
[a/Fe]. In FigureHwe showed that the degeneracy between 
[Fe/H] and [a/Fq is more mild than for age vs. [Fe/H[. 
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Nevertheless, it does exist and would have the following 
effect in our case: if the real intermediate age stellar pop¬ 
ulations of ACC-ETGs were even more [Fe/H]-poor than 
allowed by our model range, the fit would compensate by 
lowering the [a/Fe] further. We would thus end up exactly 
in the lowermost left corner of our model coverage - which 
is what we see. Our analysis therefore tentatively indicates 
that the late mass-growth of ACC-ETGs is indeed domi¬ 
nated by low-mass, low-metallicity dwarf galaxies. 

8.4. The early star formation event in SLF-ETGs 
vs. ACC-ETGs 

Another interesting feature of Figure [24] are the properties 
of the older stars. The [a/Fe]o and [Fe/Hjo are both higher 
for ACC-ETGs than for SLF-ETGs. This implies more effi¬ 
cient, early star formation for ACC-ETGs. Indeed, to reach 
such high values in both [a/Fe] and [Fe/H] a strong star- 
burst is needed, which may have been sufficient to drive out 
the remaining gas in the galaxies. Taking the speculation 
a bit further, one could imagine that ACC-ETGs would be 
different from SLF-ETGs not in the amount of material 
brought in by mergers, but rather by the efficiency of their 
early star formation event. 

Figure shows that the velocity dispersions of ACC- 
ETGs are on average higher than those of SLF-ETGs. This 
would imply a deeper potential well, which might have two 
effects: 1) increasing the efficiency of star formation, and 
2) increasing the speed needed for wind-driven gas to leave 
to leave the galaxy. Whether it can be understood through 
simulations what causes effect 1) to dominate would war¬ 
rant follow-up work. 


9. Conclusions 

We have analyzed a spectroscopic sample of 2286 early-type 
galaxies from the SDSS using state of the art stellar popula¬ 
tion models and fitting procedures. In particular the mod¬ 
els we use predict the effects of age, [Fe/H] and [a/Fe] on 
the spectra of SSPs over significant wavelength range. This 
allows us to use the code paradise to fit for the stellar 
enrichment history of the target galaxies without any prior 
assumption on any functional form or smoothness of this 
SEH. 

We conduct extensive simulations of our fitting rou¬ 
tine and find that we can reliably recover the luminosity- 
weighted average properties of galaxies. We also find that 
we can reliably recover the [a/Fe] and [Fe/H] of the 
stars older and younger than 9.5 Gyr separately. We term 
those properties [Fe/H[o and [a/Fe[o for the old stars 
and [Fe/H[i and [a/Fe]i for the intermediate age stars. 
We introduce a classification of galaxies into those where 
[Fe/H[o < [Fe/H[i - A and term those SLF-ETGs be¬ 
cause their enrichment histories are consistent with self¬ 
enrichment. AGG-ETGs on the other hand are those where 
[Fe/H[o > [Fe/H[i + A. For this paper we choose A = 0.15. 
We conclude from our study that: 

— We confirm earlier work in that the [Fe/H] and 
[a/Fe] parameters are correlated with the mass and ve¬ 
locity dispersion of ETGs. However, we find that the 
strongest relation is between [a/Fe] and age, as expected 


from basic chemical evolutiorEl This relation falls into 
two regimes, one with steep slope of 0.03 for galaxies 
with average ages above 9 Gyr and one with a shallower 
slope of 0.01 for galaxies with average age below 9 Gyr. 

— We empirically find both kinds of galaxies, i.e. SLF- 
ETGs and AGG-ETGs in our sample. For our sample, 
bona-fide SLF-ETGs make up only 8% of our sample. 
However, we find that for all galaxies independent of 
type [a/Fe] o> [a/Fe] i, as expected from basic chemical 
evolution. SLF-ETGs are on average younger and less 
a-enhanced than AGG-ETGs. 

— When studying the old stars in ETGs separately, we 
observed very high [a/Fe] ratios (above 0.3) and high 
[Fe/H] ratios (above 0.). On the other hand, and within 
our uncertainties, we find no evidence for a dependence 
of the properties of these old stellar populations on any 
other properties of their host galaxies, be it stellar mass, 
velocity dispersion or environment. 

— ETGs on average differ strongly in their enrichment his¬ 
tories from any stellar system in the Milky Way, even 
from the bulge. In particular, the vast majority of them 
does not show a ’knee’ in the plot of [Fe/H] vs. [a/Fe], 
commonly observed in local group galaxies. This implies 
that for the vast majority of ETGs the old and the inter¬ 
mediate stellar populations are not bound together by a 
common enrichment history, or in other words, the stars 
younger than 9.5 Gyrs are likely to have been accreted 
through minor mergers. Whether the accreted material 
came in the form of stars or gas (that later formed stars) 
cannot be distinguished easily. 

— The properties of the intermediate age stars in AGG- 
ETGs indicate that mass growth through late (minor) 
mergers is dominated by galaxies with low [Fe/H] and 
low [a/Fe], therefore presumably also low mass galaxies. 

— We have physically and technically motivated the choice 
of boundary between old and intermediate age stars at 
9.5 Gyr. However, the relative percentages of SLF and 
AGG-ETGs obtained in this paper do depend on that 
choice and further studies should be undertaken to de¬ 
termine the most physically and technically appropriate 
value. 

This contribution is the first one to attempt to resolve 
the stellar enrichment history of galaxies from integrated 
spectra. The method is both powerful and subject to un¬ 
certainties and caveats. We expect that further work in the 
future will allow us to confirm, extend and refine our re¬ 
sults. This will allow us to compare the SEHs obtained for 
nearby galaxies with those of much larger samples, thus pro¬ 
viding much needed comparison points and contributing to 
a united understanding of galaxy formation and evolution 
across all Hubble types. 
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Appendix A: Fitting procedure 

To derive the star formation histories of composite stellar 
population spectra we use a full spectrum fitting software 
called par adise. This software is based on algorithms de¬ 
scribed in Rix & White _|1992 1 and an earlier version was 
used previously in e.g. Walcher et al. (2006 2009 and ref¬ 
erences therein). Since then a number of new features have 
been added and we describe here the entire algorithm for 
the sake of definiteness. 


— Read in observed spectra and the corresponding stan¬ 
dard deviations. Rebin to logarithmic wavelength basis 
if original spectra are in linear. The set of observed spec¬ 
tra is considered to be in row-stacked spectrum format, 
with one axis being wavelength and the other axis allow¬ 
ing to fit any reasonable number of spectra sequentially. 


— Read in template spectra, rebin to same wavelength 
solution as object spectrum. Template spectra can be 
chosen freely, those used in this paper are described in 
Section Spectra should come with a description file 
containing information such as age, mass, luminosity, 
[Fe/H], [a/Fej. 

— Specified pieces of the spectra can be masked from the 
fit, either in restframe (e.g. emission lines) or in observed 
frame (e.g. telluric features). 

— The continuum of the object and the templates can ei¬ 
ther be left untouched, or taken out. In the second case 
the continuum can be determined through a set of Leg¬ 
endre polynomials or through a running mean with to- 
be-specified width. The continuum can be subtracted or 
the spectra can be divided by it. For continuum compu¬ 
tation reasons, a separate restframe mask can be sup¬ 
plied. 

— Use an initial template to obtain a reasonable value 
for the velocity and velocity dispersion of the object 
through a x^-biased random walk. The line-of-sight ve¬ 
locity distribution is assumed to be Gaussian. This losvd 
is applied to all template spectra. 

— Apply a non-negative l east squares algorithm (nnls, 
Lawson & Hanson|jl974|) to perform the inversion onto 
the stellar population model (see 
for details on this step). 

— Finally, repeat y^-biased random walk to yield a best 
estimate for the kinematic parameters, starting from the 
best fit values obtained before. 


Walcher et al. 2006 


The nnls routine has the feature of preferring some tem¬ 
plates over others, leading to spurious non-continuous 
star formation histories. Therefore one can choose to re¬ 
peat the stellar population inversion many times, where 
the object spectrum is noised assuming Gaussian noise 
and the variance vector, the order of the templates is 
reshuffled and a certain percentage of the template set 
spectra is not considered anymore. This is what we call 
the stellar population bootstrap. The stellar population 
bootstrap allows also to derive useful rms uncertainties 
on any of the stellar population parameters output by 
the code. In practice we find that after 100 bootstrap 
runs and using 80% of the initial template set we at the 
same time obtain smooth SFHs and useful errorbars. 
We believe this approach to be a straightforward alter¬ 


native to regularization (Ocvirk et al. 2006 Gappellari 


& Emsellem 20041. 


It is often unclear to what extent stellar population 
uncertainties can bias the velocity dispersion measure¬ 
ments. To address this, the code has the capability to 
repeat the kinematic fit after each step of the stellar 
population bootstrap. We call this the kinematics boot¬ 
strap. This feature is not used in the current paper. 
Finally the code puts out the fiducial best fit spectrum, 
the continuum it used to normalize the object spectrum, 
the associated errorbar, and residual. It also provides 
the user with two tables, one with the fiducial values of 
the fit and the other with the results of the bootstrap. 
Finally, the variance vector on the best fit derived from 
the bootstrap procedure is also provided. In this paper 
we generally use the stellar population parameters ob¬ 
tained from the bootstrap as these come with errorbars. 
Within the errorbars, these are consistent with the fidu¬ 
cial values. 
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For a review of other available spectral fitting codes, the 


f e _ 

(20111 and the website 

http: //www.sedhtting.org, 

A general note on averaging may also be useful for the 
reader: when quoting the mean [Fe/H] and [a/Fe] abun¬ 
dances, it is not trivial to determine the most appropriate 
way to do the averaging. One could either 1) average on 
the "square bracket quantities", or 2) average on actual ele¬ 
ment abundance ratios, or 3) average on actual masses of Fe 
atoms M(Fe) and of Hydrogen atoms M(H) (alpha elements 
and Fe respectively for the case of [a/Fe]). We have decided 
to stick with possibility one for an observational reason: in¬ 
terpolating linearly two stars with [Fe/H]=-1 and [Fe/H]=0 
exactly reproduces a star with [Fe/H|= -0.5 (tested using 
stellar spectra from Coelho et al. (2005)). Because we can 
thus observationally not distinguish between these last two 
situations, we shall use an averaging procedure in which 
they both correspond to the same result. In this way we 
avoid inconsistencies and degeneracies. When attempting 
to compare our results to simulations, one will have to do 
the same to their derived quantities. 
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